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^ ! Abstract 
^ '. 

^ I In the R+aR^ gravity theory, we show that if freely propagating massless particles have 

^ . an almost isotropic distribution, then the spacetime is almost Friedmann-Robertson- 

^ ' Walker (FRW). This extends the result proved recently in general relativity (a = 0), 

! which is applicable to the microwave background after photon decoupling. The higher- 

^ ■ order result is in principle applicable to a massless species that decouples in the early 
universe, such as a relic graviton background. Any future observations that show small 

. r-i '. anisotropies in such a background would imply that the geometry of the early universe 

^ . were almost FRW. 



1 Introduction 

Recently it has been shown [1] that the observed almost-isotropy of the cosmic mi- 
crowave background radiation (CMBR) implies that the universe since photon decou- 
pling is almost spatially homogeneous and isotropic. This is proof of the stability under 
perturbation of the original exact Ehlers-Geren-Sachs (EGS) theorem [2], and shows 
that almost-isotropy of the CMBR is the foundation of the assertion that the universe 
is almost a Friedmann-Robertson-Walker (FRW) spacetime. (Indeed the analysis of 
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the CMBR via the Sachs- Wolfe effect begins with the assumption that the universe is 
almost-FRW since decouphng.) 

In this paper we extend the result of [1] to a radiation-dominated universe where 
gravity obeys the higher-order theory with Lagrangian C — R + aE?. In an earlier 
paper [3] we showed that the exact EGS theorem holds also in the a 7^ theory, i.e. 
if the CMBR is exactly isotropic for a congruence of freely falling observers, then in 
R + aR^ gravity the spacetime is FRW. 

The corrections introduced by the aR^ term in the Lagrangian should be dominant 
in the very early universe. In fact the higher-order term can generate inflationary ex- 
pansion without invoking an extraneous inflaton field [4,5]. This paper investigates the 
effect this term may have on massless particles which decoupled in the early universe. A 
possible application is the relic graviton radiation, which decoupled around the Planck 
time and has since not interacted significantly with matter or radiation. In principle 
the gravitational background radiation provides a link to the spacetime geometry of the 
very early universe - although in practice any observation of this background radiation 
seems a very distant possibility. 

This radiation would inevitably contain anisotropics that carry an imprint of per- 
turbations at the decoupling time. These anisotropies would be preserved in the sub- 
sequent free propagation of the radiation. The result proved here shows that if any 
gravitational background radiation is observed in the late universe to be almost isotropic, 
and if gravity is governed by the higher-order theory soon after Planck time, then the 
spacetime geometry at that time is almost-FRW. If the higher-order corrections are ne- 
glected, the result holds for the special case of general relativity {a — 0), and represents 
an extension of [1] from a matter-dominated to a radiation-dominated universe. 

A second application, where the observational possibilities are much greater, but 
the significance of higher-order effects much less, is a (massless) neutrino background 
after decoupling. Anisotropies in this background would carry information about the 
spacetime geometry at t ^ Is. An estimate of the relative size of the higher-order 
terms at neutrino decoupling depends on the limits placed on the coupling constant a. 
In [5], a cosmological model is constructed with an aR'^ inflationary epoch, followed 
by an oscillatory phase (when re-heating takes place), which leads into the standard 
Friedmann radiation era. In this model, the inflationary constraints (e.g. on density 
perturbations) lead to the limits 

10^2 GeV^ < a-^ < 10^^ GeVl 
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According to these limits, we calculate that the relative correction to the Einstein 
Lagrangian at the start of the Friedmann era is at most 

aR 10-^^ 

This is effectively negligible, although the value is model-dependent. It is more realistic 
to take general relativity as holding at neutrino decoupling. Then our result for a = 
shows that an almost isotropic neutrino distribution implies an almost FRW spacetime 
at neutrino decoupling. 

There may be other massless species that decoupled before neutrinos, interact very 
weakly, and have yet to be detected. Finally, there is also a more mathematical motiva- 
tion for proving an almost-EGS theorem when o; 7^ 0. We can consider the higher-order 
Lagrangian as a 'perturbation' of the general relativity Lagrangian, and we show that 
the almost-EGS result is 'stable' under this kind of perturbation, having previously 
shown [3] that the exact EGS result is. 

For convenience, we summarise the theorems derived in [1-3]. Consider a congru- 
ence of freely falling observers in an expanding universe, measuring freely propagating 
radiation (massless particles). 

(i) EGS Theorem: If the radiation is exactly isotropic, then in general relativity the 
spacetime is exactly FRW [2]. 

(a) Almost-EGS Theorem: If the radiation is almost isotropic in a matter-dominated 
universe, then in general relativity the spacetime is almost FRW [1]. 
(in) Higher-order EGS Theorem: If the radiation is exactly isotropic, then in R-\-aFI? 
gravity the spacetime is exactly FRW [3]. 

In section 4 we prove: 
(iv) Higher-order Almost-EGS Theorem: If the radiation is almost isotropic in a 
radiation-dominated universe, then in i? -|- aB? gravity the spacetime is almost FRW. 

The proof of (iii) above rests on showing that the scalar curvature is spatially ho- 
mogeneous. This follows in [3] from the field equations only after detailed calculation. 
Consequently, the proof of (iv) rests on the requirement that the Ricci scalar curva- 
ture is almost spatially homogeneous. This is the main part of the generalisation of 
the general relativity theorem and is shown in sections 4.1 and 4.2. Consequently we 
can invoke theorem (iii) in our proof that the metric is almost isotropic and spatially 
homogeneous in section 4.3. In order to make the discussion as self-contained as possi- 
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ble, we briefly review in section 2 the necessary facts from the covariant formahsm for 
analysing small anisotropies in radiation [1], and in section 3 we summarise the main 
points from [3] about the higher-order equations. 

Notation: We follow [1]. The metric Qab has signature (—,+,+,+). Einstein's 
gravitational constant and the speed of light in vacuum arc 1. Round brackets on in- 
dices denote symmctrisation, square brackets anti-symmctrisation. Va is the covariant 
derivative defined by Qab- Given a four-velocity -u", the associated projection tensor is 
hab = gab + u-aUb, and the comoving time derivative and spatial gradient are 

Qa...b = u'^^cQa...b, 
^ycQa...b = K^K'...hb^VdQe...f 

for any tensor Qa...b- Given a smallness parameter e, 0[N] denotes 0{e^) and A B 
means A- 5 = [2]. 

2 Covariant analysis of radiation anisotropy 

Following [1] , we do not assume a background metric, but start from the real spacetime 
with almost isotropic radiation, and proceed to show that the real metric is close to an 
FRW metric. Of course, it follows from the higher-order EGS theorem (iii) that if the 
anisotropy vanishes, then the metric is FRW. 

We assume that the universe is radiation-dominated at the time of decoupling of 
the massless particles whose distribution is almost isotropic. The process of decoupling 
leads to anisotropies in the decoupled species, and probably also in the other massless 
or ultra-relativistic species. After decoupling, the freely propagating species preserves 
its anisotropy, while any anisotropy in the other species is rapidly removed by collisions. 

A unique physical four-velocity is defined as the unit normal field to the surfaces 
{/i = constant}, where /i is the energy density of the thcrmalised massless (or ultra- 
relativistic) species. The distribution of the decoupled species will appear as almost 
isotropic to 'observers' comoving with u"-. This field is irrotational, but in general not 
geodesic. Inspection of the proof in [1] and of the higher-order equations in [3] shows 
that it is easier to work with a geodesic but rotating four-velocity m", which may be 
chosen close to on an initial surface. Then the difference between the two vectors 
will remain provided that the spatial gradient of the energy density Jx relative to 
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is 0[1], which is proved in section 4 (see (24)). So the anisotropy measured by freely 
falhng 'observers' comoving with -u" is also small. Thus, although -u" is not uniquely 
and physically defined, the 0[1] nature of anisotropies relative to ti", and conclusions 
based on that, are covariant. 

The congruence is necessarily expanding since it must 'track' the expanding S". 
Thus satisfies 

a" = -u" = , e = u^.a > . (1) 
The rate of expansion defines a Hubble rate H and average scale factor S hy Q — 
3H — 3S/S. The kinematics of the congruence is determined by © and the shear 
and vorticity ujab- 

Now u"' defines an invariant 3+1 splitting of tensors [1,6]. In particular, for a 
massless particle four-momentum 

where E is the energy and Ca the direction of momentum, relative to the 'observers' 
above. After decoupling, the total radiation distribution function is 

ftotix'', E, e'^) = /(a:^ E) + f{x^, E, e") 

where / is a collision-dominated Planckian distribution describing the thermalised 
species, while / is the distribution function of the decoupled species. This distribution 
function may be expanded as [1,7] 

/(x^ E, e") = F(x^ E) + F„(x^ £;)e« + F„,(x^ E)e''e^ + . . . (2) 

where the covariant multipole moments Fai-.-ail^^*^, E) for L > 1 are symmetric trace- 
free tensors orthogonal to it", that provide a measure of the deviation of / from exact 
isotropy (as measured by m"). If the decoupled radiation is almost isotropic, then [1]: 

= 0[1] (L>1). (3) 

Energy integrals of the first three moments define the decoupled radiation energy 
density, energy fiux and anisotropic stress: 

fCX) 

II^Att E^FdE = 3p = O[0] , (4) 
io 

4.77- roo 

qa^^ E'FadE = 0[1] , (5) 
3 Jo 

§77- fOO 

TTab^— E^F,bdE^O[l]. (6) 
15 Jo 
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(In [1] iiR is used for /i.) We will also need the integral of the octopole moment: 

Uc = ^ / E^FabcdE = 0[1] . (7) 
oo Jo 

The total radiation energy-momentum tensor is 

Tab = (a* + fJ')UaUb + |(a* + fJ^hab + T^ab + 2li(a?6) ■ (8) 

Since the decoupled and thermalised species are non-interacting, they separately obey 
the conservation equations: 

{UT + t^Q = 0, (9) 

/i + |/xe + TTafeo;"'' +3 Vag" = 0, (10) 

ga + (c^a& + fTa6 + |e/l„5)g''+|('V„/x)+=^ = 0, (11) 

where we have used (1). 

After decoupling, the radiation obeys the equilibrium Boltzmann equation 

L{ftot) = , 

where L(/) vanishes because of detailed balancing of the collisions in the thermahsed 
component, and L{f) vanishes because the decoupled species is collision-free: 

Consequently and iTab in (8) are not dissipative quantities, but measure the deviation 
of / from isotropy. From (3), (5) and (6), it follows that 

Qa , ^VbQa - 0[1] , nab , c^^ab = ^[l] . (13) 

The Liouville equation (12) may also be covariantly decomposed into multipole mo- 
ments. The monopole moment is (10), the dipole moment is (11), and the quadropole 
moment involves ^ahc and is given in section 4 in linearised form (see (19)). 

3 Higher— order equations 

The field equations and Ricci and Bianchi identities in higher-order gravity are given 
in general 3+1 form in [3]. Quantum corrections to the gravitational Lagrangian which 
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yield the simplest higher-order field equations are of the form 



where a is a coupling constant. The field equations derived from this are 

Rab — \Rgab + 2q; [-R(-Ra6 " \Rgab) + Qab^R — R;ab — Tab , 

where □ = g'^^VaVft and Tab is the radiation energy-momentum tensor (8). Their 3+1 
splitting is 



RabU^'u^ 
RabU"'h^c 
Rabh"'ch^d 



= {l + 2aR)-^ {p + fj,) - ^R{l + aR) + 2anR + 2aR.abu''u'' 



(1 + 2aR)-^ [-Qc + 2aR.abu''h'' 
(1 + 2aR)-^ + + + aR) - 2aDR'^ hcd+ 

+7rcd + 2aR.abh''eh''d] , 



while their trace is (using (8)) 



R = GaOR 



(14) 



Together with the Ricci and Bianchi identities, these give the constraint and evolution 
equations governing a radiation spacetime in higher-order gravity. These equations are 
given in full in [3], and will be quoted when needed. 
The FRW metric is 



dr'' 



(1 -A;r2) 



where k = 0, ±1. The field equations for this metric are [5] 



2^ {S^ + k) 
3(,g^ + k) 



= -|(/i + //)-2 



a 



'S (S^ + k)' 



S 



+ R + 



2SR 



S 



{11 + li) - 2a 



' \ (^^ + kf \ 3SR 
S'']^ s 



, (15) 



(16) 



Equation (16) is the higher-order Priedmann equation for radiation and is compatible 
with (15) through the conservation equations (9), (10). In general relativity, the FRW 
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spacetime is uniquely characterised by the existence of a geodesic, expanding four- 
velocity with zero vorticity and shear. In higher-order gravity, this is also true for 
radiation, but not in general (i.e. not without further conditions) [3]. The four- 
velocity is then the unique four-velocity normal to the spatial hypersurfaces and with 
respect to which the radiation is exactly isotropic. 

4 Almost— EG S theorem in higher— order gravity 

Following the approach of [1], and given theorem (iii), we must show that if the devi- 
ations from isotropy in / are small, i.e. if (3) holds, then: (a) all covariant quantities 
that vanish in FRW spacetime are small; (b) the metric can be put into a perturbed 
FRW form. 

It is clear that any covariant quantity which is non-zero in the exactly isotropic 
(and therefore FRW) case is an O[0] quantity. However, if a covariant quantity vanishes 
in the exactly isotropic (FRW) case, it is not necessarily 0[1\. That is precisely what 
has to be proven. 

We follow the arguments of [1], explaining the deviations in the proof of the the- 
orem caused by the higher-order equations. We will only summarise the proof where 
the equations are independent of the field equations used, and include details where 
they are not. 

4.1 Almost-FRW kinematics 

The nature of the radiation kinematic quantities is independent of the field equations 
and follows from the covariant multipole decomposition of the Liouville equation (12). 
The zero and first moments are the energy and momentum conservation equations (10) 
and (11) respectively. The momentum conservation equation (11), together with (13), 
implies that 

^V„// = 0[1] , (17) 

and any of its derivatives are also at least 0[1]. The definition of leads to the 
following identity [1]: 
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By (3) and (17) this implies that 

COab = . (18) 

The quadrupole moment of (12) is the evolution equation for the anisotropic stress 
tensor TTab, which has linearised form [1] 

<b + |©7r„6 + ^fXaab + 2 {'V<ag6>} + 'Vcfab ^ , (19) 

where the angled brackets denote the spatially projected, trace-free and symmetrised 
part of the enclosed indices, i.e. 



Scd ■) 



^<ab> 

for any Sab- Using (3), (13) and (18), (19) implies 

(Tab = 0[1] , (20) 

and similarly for all of its derivatives. The linearised conservation equations are thus 

fi+lliQ + ^Vaq" ^ 0, (21) 
ga + |eg„+|(3V„/x) + 3Vft7r^ ~ 0. (22) 

Taking spatial gradients of (21) and using (4), (5) and (17), we obtain 

3V„e = 0[1] . (23) 

Equation (23) and the spatial gradient of (9) then lead to 

^Vafl = 0[1] . (24) 

Equation (23) is crucial to proving almost-homogeneity of the scalar curvature. The 
higher-order RbcU^h'^ field equation is the constraint equation [3] 

/i"^(|e;6 - abc-dh^) - v'''"^UbUJc-d = (1 + 2aR)-^[q'' - 2ah'''> RbcU''] , 

where a;" — ^rj^'^'^UbUJcd and we have used (1). This equation implies, by (5), (18), (20) 
and (23), that 

ha'RbcU" = 0[1] . (25) 



Equation (25) is insufficient to show that the spatial gradient of the scalar curvature is 
0[1]. This requires the higher-order Raychaudhuri equation [3], whose linearised form 
is 

e + + (1 + 2aR)-\jl + At) - a(l + 2aR)-^ + OR - 2i?.„b^/V] ~ , (26) 

where we have used (18) and (20). We take the spatial gradient of (26), and use (14) 
together with (17), (18), (20), (23) and (25), to obtain 

a(2e + fe^ -R)-l\ ^VcR + (1 + 2aRfVc{Q + 1©^)+ 
h/R;abdU''u'' + R.ab{hM.du' + u'^h/u'-d)] + 'Vc(// + pi) 



2a 



0. 



This implies, using (17), (18), (20), (23), (24) and (25), that 



a(|e + fe" - R) 



'VcR = 0[l]. 



Since they are non-zero in the exactly isotropic (and thus FRW) case, R, © and © are 
O[0]. Hence we deduce that 

^VaR = 0[1] and thus A<ed> = 0[1] . (27) 

The definition of the magnetic part of the Weyl tensor Hab in terms of the kinematic 
quantities (18) and (20) is independent of the field equations so that, as in [1] 

Hab = 0[1]. (28) 

However, to prove that the electric part Eab of the Weyl tensor is 0[1] requires the 
higher-order shear propagation equation [3]: 

- i(l + 2ai?)-V„6 - a{l + 2aR)-^ R,<cd> = , (29) 
on using (1). Then, by (18), (20), (27) and (29) 

Eab = 0[l] . (30) 

This establishes that the kinematic quantities for the radiation and geometry are 
almost-FRW, i.e. all the covariant quantities that vanish in the exactly isotropic (FRW) 
case are at least 0[1\. The higher-order terms in the field equations thus have no O[0] 
effect on any of the kinematic quantities. 
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4.2 Almost— FRW dynamics 

It follows from section 4.1 that the O[0] equations governing the dynamics are just those 
of an FRW universe. Thus by (26) and the conservation equation (21), the evolution of 
the scale factor S is to O[0] that of a higher-order radiation FRW universe. The first 
integral of (26) is the almost-homogeneous analogue of the higher-order Friedmann 
equation (16) 

Then the dynamical equations differ from those of the exactly isotropic (FRW) metric 
at 0[1] only. Because the kinematic quantities are just those of a perturbed FRW 
universe, it is now possible to linearise the covariant equations about the FRW values 

in the usual manner [8], and to consider the evolution of density inhomogeneities. The 
background metric will obey the FRW equations, and consequently we can use the 
usual covariant FRW perturbation analysis. 

4.3 Almost-FRW metric 

By returning to the frame u"- described in section 2, we now show that the metric 
may be given in an almost-FRW form. In [1], the choice of is based on surfaces of 
constant matter density, since for applications to the CMBR, the universe is matter- 
dominated. In our case, is based on the surfaces of constant thermalised radiation 
density {ua oc p,,a), and 

= , a V . 

(All kinematic quantities of the -u" congruence are denoted by a hat.) 

We are able to choose the geodesic to be close to because of (24), which 
ensures that the angle between and it" will be 0[1]. Clearly the energy flux and 
anisotropic stress relative to remain 0[1\. Changes to all 0[1] quantities will be 
0[2], but changes to O[0] quantities will be 0[1], and will introduce effective energy 
fluxes with respect to S". 

The congruence will be normal to cosmic time surfaces (coinciding with surfaces 
of constant p) in the perturbed model, and the spacetime metric will be [1]: 

ds^ = -A\x'')dt^ + Kpdx'^dx^ , (32) 



g2 




, k^O[l\. (31) 
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where = A ^(5"o , hab = gab+UaUa, and Greek indices run from 1 to 3. The kinematic 
quantities for fi" are [1] 

a« = = -h''\\ogA),b , cD„(. = , = lAhab,o • (33) 

If we express the non-coincidence of the flows of u"' and by [1] 

S» = li" - T/" + , UaV'^O , K^" = 0[1], (34) 

then the kinematic and dynamic quantities relative to are given in terms of those 
relative to by 

e-e + 0[l], aab^aab + 0[2], a« = 0[l] 
= + g« = g» + |//y'^ + 0[2] , nab ^ TTab + 0[2] , (35) 

where p, is the energy density of the decoupled species as measured by 'observers' 
comoving with the thermalised species. Thus the kinematic results (17-30) hold in 
this frame, and in particular 

V„/i = 0[l] , V„e = 0[l] . (36) 
Since the expansions are equal (up to 0[1]), the scale factor S{t,x"), defined by 

A-\logS),o=le : 

will satisfy the higher-order Friedmann equation (31). From (36) and this definition 
of S we can write the spatial metric in (32) as 

ha^{t, X^)dx''dx^ = S\t, X"')fa^{t, X^)dx''dx^ , (37) 

where most of the time variation is in S, and from the result [1] 

^ab = ^AS'^fabfi = 0[1] , 

most of the spatial coordinate dependence is in fab- Equations (33) and (35) imply 
that the function A(t, x") is only weakly dependent on the spatial coordinates, and we 
can consequently set AitjX"') = 1 + 0[1] by rescaling the time coordinate. 
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It remains only to show that the embedding of the 3-surfaces is almost-isotropic. 
The higher-order Gauss-Codazzi equations [3] for the surfaces {jl = const} give the 
trace-free part of the Ricci tensor of these surfaces as 



-1 



^Rab — \^R)hab — hg^hb'^ 

(1 + 2aR)-^ l7f„b + 2a{haV - \habh"')R;cd 



(38) 



The scalar curvature is independent of the four-velocity, and consequently (27) gives 



(hSb" - \habh'^) R-cd = 0[l] . 



(39) 



Given (39) and the almost-FRW conditions on the kinematic quantities (35), equation 
(38) implies 



Rab-U'R)Kb = 0[l\ , 



(40) 



where, from [3] 



R = -|e2 + (1 + 2aR)-^ [2/2 + aR^ + AahP-^R.J + 0[1]. 



(41) 



Because R, and /2 are almost spatially homogeneous, by (36) and (41), so is ^R. 
Then the 3-spaces are spaces of almost-constant curvature: 



^ Robed — 



hachbd — hadhbc 



+ o[i] . 



Thus, as in [1], we have recovered (to O[0]) all of the standard relations governing an 
FRW universe, and in particular, we have shown that there exists an almost-FRW 
metric in higher-order gravity, given in almost-comoving coordinates. 



We have shown that in a radiation-dominated universe, an almost-isotropic back- 
ground radiation implies an almost-FRW spacetime geometry in the higher-order the- 
ory of gravity (including general relativity as a special case). In effect, together with 
the higher-order exact EGS theorem of [3] , we have now shown that the EGS theorem 
is stable under perturbations of both the Lagrangian and the background metric. 
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